FINITE VOLUME SCHEMES OF ANY ORDER 
ON RECTANGULAR MESHES 
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Abstract. In this paper, we analyze vertex-centered finite volume method (FVM) of any order 
for elliptic equations on rectangular meshes. The novelty is a unified proof of the inf-sup condition, 
based on which, we show that the FVM approximation converges to the exact solution with the 
optimal rate in the energy norm. Furthermore, we discuss superconvergence property of the FVM 
solution. With the help of this superconvergence result, we find that the FVM solution also converges 
to the exact solution with the optimal rate in the L 2 -norm. Finally, we validate our theory with 
several numerical experiments. 
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1. Introduction. During the past several decades, the finite volume method 
(FVM) has attracted much attention. We refer to PJ-0, P3]-[l2], [22]-[2S]' [2Qj[35] for 
an incomplete list of references. Due to its local conservation of numerical fluxes and 
other advantages, the FVM is very popular in scientific and engineering computations, 
especially in computational fluid dynamics, see, e.g., [13l [TTl [18] and [22" ] -[26 ] . 

Comparing to its wide applications, the mathematical theory of FVM (cf., [3] 
[15l [19l [20]) has not been fully developed, at least, not as satisfactory as that for the 
finite element method (FEM). In fact, since the FV schemes depend heavily on the 
underlying meshes, the error analysis in the literature was often done case-by-case. 
For instance, the linear FV scheme can be regarded as a small perturbation of its 
corresponding linear FE scheme, whose convergence properties have been well studied, 
see e.g., [2] [4] [Ml [16]. On the other hand, high-order FV schemes are substantially 
different from their corresponding FE schemes, therefore only a few special high-order 
schemes have been studied, see [5J [7] El EU [25l [27l [30] . So far, we have not seen analysis 
for FV schemes of an arbitrary order. 

In this paper, we provide a unified analysis for vertex-centered FV schemes of any 
order on rectangular meshes. We construct our FV schemes under the framework of 
the Petrov-Galerkin method by letting the trial space be the Lagrange finite element 
space with the interpolation points being the Lobatto points and by constructing 
control volumes with the Gauss points in a rectangular element. 

It is known that the proof of the stability is a challenging task in the error analysis 
of FV schemes. Earlier works (see, e.g. [21] (20] [30] [8]) utilized element stiffness 
matrix analysis for this task. The element stiffness matrix analysis often requires 
to calculate all eigenvalues of an element stiffness matrix and thus is difficult to be 
generalized to schemes of any order. Our new approach in this paper for proving the 
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stability (or in general the inf-sup condition) is different from the element stiffness 
matrix analysis. A novel and non-traditional global mapping from the trial space 
to the test space is introduced. This mapping avoids calculating eigenvalues of an 
element stiffness matrix and makes the establishment of the global inf-sup condition 
for any order possible. An interesting feature is that when the coefficient a in (I2.1|) is a 
piecewise constant function with respect to the underlying mesh, the inf-sup condition 
is uniformly valid for any mesh size h, i.e., there is no requirement "for sufficiently 
small h" . In particular, for the Poisson equation, the inf-sup condition is uniformly 
valid for any mesh size h. Once the inf-sup property has been established, the error 
analysis in the energy norm is then a routine work. 

Another feature of this work is the superconvergence analysis. We prove that the 
FV solution u-p is super-close to the Lobatto interpolant Uj of the exact solution, 
namely, \u-p — uj\i converges one order higher than the optimal rate. The result simu- 
lates the counterpart result in the FEM. A by-product of this superconvergence result 
is the optimal L 2 error estimate. Conventionally, the L 2 error estimate is accom- 
plished by the duality argument or the so-called Aubin-Nitsche trick. Unfortunately, 
this technique is very difficult to be used in our case for higher-order FVM. The 
adoption of the superconvergence analysis avoids this difficulty. 

We organize the rest of the paper as follows. In Section 2 we present FV schemes 
of any order for elliptic equations on rectangular meshes. In Section 3 we provide 
convergence analysis and establish the optimal convergence rate in both H 1 and L 2 
norms. The superconvergence property of the FVM solution has also been studied 
in this section. Next, numerical examples are provided in Section 4 to confirm our 
theory. And lastly, some concluding remarks are given in Section 5. 

In the rest of this paper, "A < B" means that A can be bounded by B multiplied 
by a constant which is independent of the parameters which A and B may depend 
on. "A ~ B" means "A < B" and "B < A" . 

2. FVM Schemes of Any Order. In this section, we present finite volume 
schemes of any order to solve the following second-order elliptic boundary value prob- 
lem 



where £1 = [a, b] x [c, d] is a rectangle, V = dCl, a £ L°° and it is bounded from below: 
There exists a constant cto > such that a(x) > «o for almost all x £ f2, and / is a 
real- valued function defined on f2. 

We present our finite volume schemes in the framework of Petrov-Galerkin method. 
We first construct the primal partition V and the trial space. Let a = xq < x\ < . . . < 
x m = b, c = yo < yi < . . . < y n = b. For a positive integer k, let = {1, . . . , k} 
and Z° = {0, 1, ... , k}. For all i € Z m ,j G Z„, let hf — Xi — Xi-i, frj = y.j — yj-i 
and h = max{max(/if ,hj)\(i,j) £ 1 m .n\ where Z m .„ = Z m x Z n . We denote the 
associated partition of by 



where t^j = [xi—i,xi\ x [yj-i,yj]. We choose the trial space as the standard FEM 
space defined by 



V • (ayw) =/ in 

u =0 on r, 



(2.1) 
(2.2) 



V = {n-j\{i,j) £ Z m) „} 



UJp = {v £ C(fl) : v\ T £ Q r ,Vr £ V,v\on = 0}, 
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where Q r is the set of all bi-polynomials of degree no more than r. Obviously, 
dim Uj, = (mr — l)(nr — 1). 

We next describe the dual partition and the test space. Let G\, . . . , G r be r Gauss 
points, i.e., zeros of the Legendre polynomial of rth degree, on the interval [—1, 1]. 
For any given G Z m , n , let 

9i,k = l&i + Xi-i + K G k ), g v u = \{Vi + Vj-i + h^Gi)- Vfc, I e z r 
be the Gauss points in the interval [xj_i,Xj] and [yj-i, yj], respectively. Let 

T if i G Z m , „y _ ( r if j e Z n , 
1 if i = 0, T i ~ \ 1 if j = 0, 

and 

5(3,1 = a i 9o,l = C ' 9m,r+l = 9n,r+l = 

and for all i G Z m _i, j G Z n _i 

With these Gauss points, we construct a dual partition 

P'=Kfc ; ;,i|(^,.7>0e£o}, 
where Zo = Z^ X Zf| x Z° x Z^h and the control volume 

The test space Vp/ consists of the piecewise constants with respect to the partition 
V 1 which vanishes on the boundary control volumes. In other words, 

Vp> = Span {ipi,k;j,l | (h k; j, I) G Z x } 

where ipi t k;j,l = Xt' fc . , is the characteristic function on the control volume t^j.. -;, 
Zi = Z m x Z r * x Z n x Z„h and 

r if i e Z m _i, r v _ [ r ^ j G Z n _i, 
r-1 if i = m, J ' 1 ?* — 1 if j = n. 

Note that dimVp' = (mr — 1) X (rtr — 1) = dimt/£,. 

We are now ready to present our finite volume schemes. The finite volume solution 
of (|2.1[) and (|2.2p is a function «p G which satisfies the following conservation law 

a n V ds = [ fdxdy (2.3) 
an ; 



on each control volume t- k .j ; , (i,k; j,l) G Z\ , where n is the unit outward normal on 
the boundary curve dr' ik .- ; . Let wj>j G Vpt, w-pi can be written as 

(i,k;j,l)eZi 
3 



where the coefficients Wf ^.j i are constants. Multiplying (I2.3[) with Wi^-j.i and then 
summing up for all i, k, j, I, we obtain 

- V] w l . k - j .i [ ot^-ds = [ fwp,dxdy. 

Defining the FVM bilinear form for all v G Hq(Q.),w-p' G Vv as 

f du 



Or' 



the finite volume method for solving equations (12.11) and (|2.2j) reads as : Find up g Uj, 
such that 

ap(up,wp>) = (f,w-pi), \/w-p> G Vp>. (2-5) 
Noticing that Zi C 2p, a function wpi G Vp» can also be written as 

WV' = X! w i,k;j,l^i,k;j,h 
(i,k;j,l)eZ 

if we let 

m,k-,j,i = 0; V(i, fc; j, Z) G Z \ Zi . 

For all (i,k;j,l) e Z2 = Z m x Z r x Z„ x Z r , we define the jumps of w-p' in the 
cc-direction and y-direction as 

[wi,k;j,l} X = Wi.fcjj',! _ w i,k-l;j,U [ w i,k;j,l] V = M i,fc;j',* ~ Wi,k;j,l-1, 

respectively. 

A simple calculation yields that 



ap(u,u)p/)= 2^ [Wi,k;j,i] X I a {9ik-,v)-Q^{9i,k,y) d V 

^ rg?,k+i Q v 

22 [ w i,k;j,t] v a ( x '9li)-Q-(x,gli)dx, (2.6) 



(i,k;j,l)eZ 4 g i 

where Z3 = Z m x Z r x Z ra x Z r y and Z4 = Z m x Z r x x Z„ x Z r . 



3. Error Analysis. The error analysis of FVM can also be done under the 
framework of Petrov-Galcrkin methods, see [2J, [2D], and [3D]- Following this ap- 
proach, our FVMs error analysis requires the study of the continuity (boundedness) 
and inf-sup property of the FVM bilinear form. 

3.1. Continuity. Let &pi be the set of interior edges of the dual partition V' . 
Then for all v G H&(Q), w v > G Vp>, 

a v (v,Wp>) = S~] [w v ,] [ a^-ds, (3.1) 
Ee£ v , Je dn 
1 



where [w-p>] = w-pi | T2 — w-p> \ T1 across the common edge E — t\ fl r 2 of two rectangles 
Ti, r 2 G T 5 ' and n denotes the normal vector on E pointing from t± to t 2 . 

To study the continuity of a(-, •) , we define a semi-norm in the test space Vp> for 
all w-pi G Vp/ by 



" = E ^ / [wvfds 
\EeS-p, Je 

where He is the diameter of an edge E, and a semi-norm in the so-called broken H 2 
space 

H%,(Sl) = C(Q) : «| T G # 2 ,Vr 6 

for all v G H^{VL) by 

\T£V / 

where h T is the diameter of r. Note that the mesh-dependent semi-norm | • | has 
been used in the discontinuous Galerkin method (cf., [T]) and was introduced first 
into the FVM in [50] . 

Theorem 3.1. The finite volume bilinear form a-p(v) is variationally exact: 

a-p(u,wp>) — (f,w-p>) V w-p' G Vp> (3.2) 

and continuous : for all v G Hq(Q) (~l Htp(Q,), w-p> G Vp/, 

|a-p(v, wp<)| < Mlvj-pluip/ |-p/ (3-3) 

where the constant M > depends only on a and r. 

Proof. First, Q3.2[) follows by multiplying (|2.ip with an arbitrary function G 
Vp' and then using Green's formula in each control volume r G V . 

Next we prove (13. 3|) . By the Cauchy-Schwartz inequality, for all v G Hq(CI), wv G 
Vp', we have 



a-p{v,w v <) < ||a||oo|u>7"|7" f E ^ B y (f^) 



2 



By the trace inequality and the shape regularity of V, 

^tlE J (J^J ds^j < |«|l, Tl UT2 + hE\v\l, Tl UT 2 

where ri, t 2 G V' are two control volumes sharing the common edge E. Therefore, 

a V {v,W V .) < \w-p,\ V ,\ ^ l^ll.nU-ra + /l ll U ll,riUT2 
\Ee£ P , 

<\W V ,\ V , [EKr + ^rMl,r 
\t£V / 

Then the boundedness (|3.3[) is proved. □ 
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3.2. Inf-sup condition. This subsection is the core of the paper. The analysis 
here is technical, and yet, it is new and no-traditional. A key step is the introduction 
of a global projection (I3.4[) based on the idea of the Gauss quadrature. We begin with 
some definitions and notations. We denote by Aj,j e Z r the weights of the Gauss 
quadrature 

r 

Q r (F) =Y J A 3 F{G j ) 

for computing the integral 

1(F) = J F(x)dx. 

It is well-known that Q r (F) = 1(F) for all F <E P 2r _i(-1, !)• For any given (i,j) e 

be the Gauss weights in the intervals [xi-\,Xi] and [yj-i, yj], respectively. On the 
other hand, [wi^k-j,i] x and [wi t k-j,i\ v are well defined for all (i, k;j, I) G Z2 and we have 

[[Wi,k;jA X ] V = [[wi,k;j.l] V ] X = W i,k;jJ + Wi,k-l,j ,1-1 ~ Wi,k-X;j,l ~ Wi,k;j,l-1- 

We denote 



lw]i,k;j,l = W i,k;j,l + W i: k-l;j,l-l ~ W it k;j,l-1 , V(«, fc? j, I) £ Z 2 . 

We define the mapping from the the trial space Uy> to the test space Vp> as : 

Uv-p = ^2 (Rv-p)i f k;j,ii)i,k\j,l e V-p',Vp S C/p, (3.4) 
(i,k,j,i)eZi 

where the coefficients (J\vp)i t k-j,i arc determined by the constraints 

IILvtUm = A;,Al^: !f ;,.,jlu :<;.k:j.h e Z 1 . (3.5) 

Remark 3.2. We can not suppose a priori that (13. 5[) holds for all (i, k;j, I) S Z2, 
because j^Zi = mnr 2 which is greater than dim V-p> — (mr — l)(nr — 1). However, 
^Z\ = (mr — l)(nr — 1), so we constraint (|3.5|) only for the indices in Z\. We present 
the Gauss points corresponding the index sets Z\ and Zi in Figure 11 (r — 2). In 
this figure, the gauss points corresponding to Z\ are depicted with ' and the Gauss 
points corresponding to Zi \ Z\ are depicted with heavy '* '. 

We next explain how to determine Uv-p by (|3.5p . For any given v-p, Uv-p = on 
the boundary control volumes, namely, (^ivp)i^k-j.i — 0,V(i,fc; j,l) £ Zq\ Z\. In the 
first step, let (i, k;j, I) — (1, 1; 1, 1) in (|3.5[) to obtain 

(n« 7 ,) 1 , lilll = Af >1 A* 1 ^(^ il>fl » 1 ). 



Once (nv-p) lil;ljl is obtained, we then let (j, I) S Z„ x Z rj \ {(1, 1)} in (I3.5[) to obtain 

(ILv-ph^j = : • Ai t g[, . g^,). 

where we again used the fact Tlv-p = on the boundary volumes. That is, (Ilv-p) 1 1 ; j t \ 
can be successively calculated for all (j, I) G Z„ x Z r . . In the same way, we can 
calculate (Hvp)i, k-i.i for all (i, k) £ Z m x Z ri successively. In the second step, we use 
(|3.5p to compute all (IIwp)i^:j,i and all (IIu-p)^;!^ for all (j, I) s Z„xZ ri \{(l, 1), 1, 2} 
and all (i, k) e Z m x Z ri \{(1, 1), 1, 2}. In the p-th step, p = 2, . . . , min(mr — 1, nr — 1), 
we use (|3~5j) to compute all (IIvp)^^..,-,; and all (n.v-p)i,k;i„,k p for all (j, Z) ee Z„ x 
Z rj \{(1, 1), 1,2,... (»„, fc p - 1)} and all (i, fc) e Z m x Z rj \{(1, 1), 1,2,..., (i p , j„ - 1}. 
Finally, we obtain (Jk&p)i t k;j,i for ah (i,k;j,l) 6 Si. 

Next, we show that (|3.5[) holds for all (i, fc; j, i) £ Z2. In fact, since vp — on the 
boundary dft, 

then for all (i,k) £ Z m]7 ., 

}2„,_ r d 



In other words, 



O^V-p x y \ \^ ty 9^ V-p x y 



Consequently, 



' 3 



= - 2 [[(n^k fe;jy ]T 

' 3 

= -[nu P ]f ife .„ r = [nupj i)fc . nir . 

By the same reasoning, for all (j, I) € Z„ iT ., 

Lnu-pJ m , r; j,i = Af llr A v jl -^-^(g^ nr ,g^ l ). 

Namely, (|3.5p holds for all (i, k;j, I) € 

Lemma 3.3. //"P is shape regular, then for any vp € Up,, 

|Hup|7" < \vp\i, (3.6) 

where the hidden constant depends only on r. 
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Proof. By the definition of the semi-norm | • |-p<, we have 

\uv v \ 2 v ,= e mwhk^rf + E ([(n^k^n 2 . (3.7) 

{i,k;j,l)eZ 3 (i,k;j,l)eZ 4 

We next estimate \(Ilv-p)i,k;j,iY ' ■ For all given (i,k;j) £ Z m xZ r x Z„, the function 

dxdy 



-^(gfki ') is a polynomial of degree r — 1 in the interval [yj-i, yj], therefore, 



zez„ y ■'s/j-i 



dxdy 



Noticing that 



we have 



dv-p 
dx 



dv-p dv v 



Consequently, for all (i, fc; j, /) £ 2 3 

(j',i')eZjxz rj . , 



.4 



E.J, (PVV x y 



i.fc ^ dxdy 1 



(3.8) 



where r^j = r, if j' < j and = Z, if j' = j 



We first estimate Since ^(fl? fc) -) £ Pr-i^-i, by the i 



verse inequality, there holds 



<9x 



dx 



Loo(.[Vj-l,Vi]) 

where the hidden constant depends only on r. Then 



L 2(lvj-i-,yj}) 



dV V , X N 



< 

L z([Vj-l,Vj]) 

We now estimate the second term of (|3. 81) . By the Cauchy-Schwartz inequality, 

2 



(3.9) 



I'EZr v 

=r ^r 1 (^ ( ^ ,2/) ) dy - 



2 / 3> 



Again since far(gfki ') ^r-iiVj-itVj), we have the inverse inequality 
d 2 v-p 



dydx 



L2 {Vj-l,Vj) 



dx {9t ^ > 



L 2 (yj-i,yj) 



where the hidden constant again depends only on r. Consequently, 



(3.10) 



Substituting (|3.9|l and (13. 10|t into (|3.8I) and noticing the fact Af k h 1 < 1, we obtain 
that for all (i, k; j, I) S Z3 

(\n»p]Zk-j t i) 2 Z (^G&>w)) d y- 



Consequently, 



( dvv ( 



E (Wvvhk-jA*) 2 < E ^f, fe / l^fo ; ''//■ 



Since the function J fi d (^§^-{ m , y))* dy 6 P2, — 2(2^-1! Xi) for all i G Z m , we obtain 



E ([(H^k^n 2 < £ £ (j£{x, 2/)) ' 



(i,ft;i,i)e^3 
By the same arguments 



(3.11) 



E (pupWl'f^^^^y)) &dy. (3.12) 



{i,fc;j,/)G24 

Therefore, by (|3~7|) 

\nv v \ 2 v , < [ ( 

J a J c 

from which the inequality (|3.6I) follows. □ 

Lemma 3.4. J/ a is piecewise constant with respect to V , then 
a-p(v-p,Hv-p) > ao\v-p\ 1: Wv-p S U^,. 



-^r (a; ' y) ) + (^ (a; ' y) j dxdy ' 



(3.13) 



Proof. We define two functions for all (a;, y) £ f2 by 

v 1 (x,y)=[ a(x,y')^-(x,y')dy / , v 2 (x,y) = [ a(x' , y)^-(x', y)dx' 



dx J a dy 

A straightforward calculation yields that 

a v (v v ,Uv v ) = - E l Uv v\t,k- 3 d (v 1 {gl k ,gl l ) + v 2 {gl k ,gl l )) 



(i,k;jA)£Z 2 



Noticing that (13.51) holds for all (i,k;j,l) E Z 2 , we obtain 

a-p (v-p , TLv-p ) = I\ + I2 , 



where 



and 



We next estimate 7i. To this end, for all (i, fc; j) £ Z m x Z r x Z n , we denote by 
Err 3 {9l k ) = £ ^(l$,k,vVUi*V)dV ~ E A l^^ k , {9^, 9^) 

the error of the Gauss quadrature of the function (fff fc; y) yl {9i k^ ') m ^ ne interval 
[j/j _ 1 , j/j ] . Moreover , let 

#es= E ^a^'(^fe). 

(i,fc;i)ez m xZ r xZ n 



With this notation, 



h = - E A h f jr^(9t,k,y)v 1 {9tk,y)dy + Res 

E / a (9i,k>y) (^(9i,k,y)) dy + Res. 



(i,fe)ez m x 

Since a is a constant in each t eV, 
then 



/1 = J J a(x,y) (-^ri x , y)) dydx + Res. (3.14) 

We next estimate Res. Using the fact that a is piecewise constant with respect to V, 
we have 

^C^fc.-y^fc.O G P2r(b/i-l,!/i]),Vj G Z„. 
Then for all y G [yj-i^j], 

Q(2r) / Q2 y \ r 

I dxdy^' k ' y ^ V ^ fc '^ ) = a ^ fc >2/)(2r-)! x x > 0. 
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where a^-j is the leading coefficient of the polynomial ^§£-(gfk, y) in [Vj-iiVj]- Con- 
sequently, 



and thus 



In summary, 



Errj(gf k ) > 0,V(i,fc) G Z m x Z r x Z„, 



i?es > 0. 



/i > «o / 6 /" f^O^ dydx. (3.15) 



By the same arguments, 



h>a J J (^j-(x,y)J dxdy. (3.16) 

Combining (|3.15|) and ()3.16| ). the inequality (|3.13|) follows. □ 

Summarizing the above two lemmas, we obtain the following inf-sup property. 
Theorem 3.5. Let V be a shape regular and quasi-uniform partition with the 
meshsize h. If the coefficient a is piecewise constant, then the inf-sup property 

. a v (v v ,W V >) 

ml sup | — ; — > 1 (3.17) 

■»T&T,v3 V ,^V v , \vv\\\wv'\v 

holds for all h. If a is piecewise continuous, then (j3. 1 7[) holds when the meshsize h 
is sufficiently small. 

Proof. When a is piecewise constant, by (|3.13[) and (I3.6[) . for any v-p G UJp, 

a v (v v ,w-p>) a-p(v v ,Uv-p) a \v v \j 
sup — ; > — — — . > — — — > Dp i. 

w v ,ev v , \wv\v \Uvv\v \Uvv\v 

The inf-sup condition (I3.17|) is proved. 

When a is only piecewise continuous, let 

\T\ Jr 

and denote the piecewise modulus of continuity of a by 

m-p(a,/i) = sup{|a(xi) — a(x 2 )| : |xi — x 2 | < h, V xi,x 2 G t,V t G V}. 

The fact that a is piecewise constant implies that m-p(a,h) converges to when h 
goes to 0. Since a is piecewise constant, by Lemma [33 



d-p(v-p,Uv v ) := y~] [Ilv-p] / a—ds > ao\v v \\. 
Jri Je on 



On the other hand, by the same arguments in Theorem 13. 1[ we have 

\a-p(v-p, Ilv-p) — dp(vp, Hvp)\ < m-p(a, h)\v-p\1, 
Then when h is sufficiently small, 

I 1 2 I 1 2 

ap(vp,Iiv v ) > (a - cm v (a, hj) \vp\ in > — \Vv\ 1M - 



The inf-sup condition (I3.17|) is proved. □ 
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3.3. H 1 error estimate. We begin with some preparations. First, by the inverse 
inequality, 

\vr\ v ~ \vv\ v Vvt> g U v . 
With this equivalence, we deduce from the inf-sup condition (|3.17p that 

. a v (vp,w v >) 

ml sup : — j : — > c , (3.18) 

Vv^Ulp Wv , e y v , \V-p\-p\W-pi\-pl 

where the constant cq = Co(r) depends only on r. 

Secondly, we denote by Li,...,L r -\, the zeros of L' r (x), where L r (x) is the 
Legendre polynomial of order r in the interval [—1,1]. Moreover, we denote Lq = 
—l,L r = 1. The family of the points Lj,j = 0, . . . ,r are called the Lobatto points of 
degree r. 

The Lobatto points on the rectangle r^j, G Z m n are defined as 

k,k;j,l = + x i-i + ^f-^fe)j \{V] + Vj-l + hV ] L i)^j ' e z n-r- 

Let Ui G f/p be the interpolation of u that satisfies 

ui(h,k-,j,i) = u(k,k-,j,i),(i,k;j,l) G Z . 

Note that this kind of interpolation has been used in the literature for superconver- 
gence analysis, see, e.g., [32l [33] for the one-dimensional situation. 

We are now ready to present our main result. 

Theorem 3.6. Let u be the solution of (|2.1[) and (|2.2[) . u v the solution of (|2.5I) . 
Then for sufficiently small h, 

M 

\u — u-p\„< — inf \u — v-p (3.19) 

' Co v v eu-p 1 1 ' 

Consequently, if u G if r+1 (fi), 

< h r \u\ r+1 , (3.20) 

where the hidden constant is independent of the mesh size h. 

Proof. By (|3.2[) . (|3.3|) and the inf-sup condition (|3.18|) . for all v-p G U-p, there 
holds 

i ii ii \ / M \\ i 
\u — u-p\ v < \u - vp\ v + \ vp — up\ v < (1 H j |u - 

Using a technique in Xu and Zikatanov ([H]), the constant 1 + ^f- in the above 
inequality can be reduced to That is, (|3TT9"|) holds. 

We conclude from the definition of | • |p and (13. 19)) that 

In — up\\ < lit — Upl-p < inf |u — v-p\v- 



Note that 



inf | it — vp\p < \u — ui\i + h\u — U7I2, 



where uj is the Lagrange interpolation of u at the Lobatto points in the trial space 
Uj, . By the standard approximation theory, we obtain the estimate (|3.20p . □ 
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3.4. Superconvergence and L 2 error estimates. We first present a super- 
convergence result and then use it to establish our L 2 error estimate. 

THEOREM 3.7. Assume that u 6 H^(Q) n H r+2 (fl) is the solution of flUD - flOD , 
and up is the solution of the FV scheme (|2.5[) . Then for all uipi G Vp>, 

\ap(u — ui, w-p/)\ < /i r+1 |w| r+ 2 i -p|u'p' |-p', (3-21) 
i 

where \u\ r +2,v = (J2 T ev \ u \r+2,r) 2 ■ Consequently, 

\ui -U V \i < h r+1 \u\ r+2 ,v- (3.22) 



Proof. We can derive the following inequality by the standard superconvergence 
argument, see, e.g., [34] . 



d(u — uj) 



dx 



d(u — uj) 



dy 



< h r \u\ r+2 ,i.f, k ., n 



where n.k.j.i = \9i,k-n 9i, k +i\ x [fl^-uS^+il- rt follows from (EHJ) that 



(3.23) 



\av{u-ui,W' P ')\<\\oi\\ 00 \'Wv'\v l \ ^ h E 



d(u — uj) 
dn 



ds 



< h r+1 \wp,\ v ,\u\ 



r+2,V, 



where in the last step we have used (|3.23j) and the fact |u| r +2.i,n ^ h£ l u k+2,r i and 
We next show (|3.22[) . By the inf-sup condition (|3.17|) . 



\ui - u-p | i 



< 



sup 



a-piu-p — Uz, w-p') 



\Wp> \pi 



sup 

w v i ev v , 



ap(u — uj, wpi) 



\wv\v 



Combining the above inequality with (|3.21j) . we derive (I3.22[) . □ 

As a direct consequence of the superconvergence result (|3.22[) . we have the fol- 
lowing L 2 estimate. 

Theorem 3.8. Assume that u £ H^(fl) n H r+2 (Q) is the solution of ([2~T ]l -([2~2 ]l . 
and up is the solution of the FV scheme (|2.5[) , then there holds 

(3.24) 



\\u-Up\\ Q <h r+1 \\ U \\ r+2 . 



Proof. By the triangle inequality, 

\\u - up\\ < \\u - ui\\ + \\u-p - u/||o 

where ui is the interpolation of u given in the previous subsection. 
Since ui = up = on dfl, we have by the Poincare inequality that 

\\up - ui\\o < \u v - < h r+1 \u\ r+2 . 

Moreover, 

h-"/||0<^ +1 h||r+l<^ + 1 || U ||,-+2. 

Then we have (pT2"4]) . □ 

Remark 3.9. In the above L 2 error estimate, we do not need to use the so-called 
Aubin-Nitsche techniques. 
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4. Numerical Results. In this section, we present a numerical example to 
validate the theoretical results proved in previous sections. We consider (I2.1[) - (|2.2I) 
with a = 1 and H, = [0, l] 2 . We choose the right-hand side function 

f(x, y) = 2n 2 sin7rx sinny, (x, y) E [0, l] 2 

which allows the exact solution 

u(x,y) = sinnx simry, (x, y) E [0, l] 2 . 

We use FV schemes (|2.5|) with r — 2,3,4,5 to compute FVM approximate solutions 
of u. The partition Vj, j — 1,...,6, are obtained by uniformly refining the unite 
square [0, l] 2 . For simplicity, we write Uj, instead of u-p., as the finite volume solution 
u Vj E . 

The numerical results are demonstrated in Figures 1, 2, 3, and 4, respectively. In 
these figures, we depict nj r by the solid curve with '□' and depict n - ' r+1 -* by the 
dash line, where rij is the number of subintervals of the underlying grid Vj- We depict 
\u — Ujlii 1 by the solid curve with V, \\u — Uj\\i,2 by the solid curve with 'A', and 
\ui — Ujlu 1 by the solid curve with '*'. We observe that \u — Uj\jji decays with the 
convergence rate n~ r which supports our theory (|3.20[) . We also observe that both 

||u — MjUl 2 an( i \ u i ~ u j\h 1 decay with rij These facts support our I? estimate 

(I3.24[) and superconvergene result (13.221) . 

5. Conclusions and Future Works. The design and analysis of high-order 
FV schemes are challenging tasks. This paper is the second one in its series that 
attempts to set up a mathematical foundation for a family of high order FV schemes. 
In a previous work ([6]), we studied convergence and superconvergence properties 
of FV schemes of any order for the one-dimensional elliptic equations. The higher 
dimensional case is fundamentally different from, and much more complicated than the 
one dimensional case. In this article, we only report our results for rectangular meshes. 
For the family of higher-order FV schemes discussed in this paper, we obtained almost 
the same basic theoretical results as for the counterpart higher-order FEM. The results 
for more general meshes and more general equations will be reported in a forthcoming 
paper [31]. 

Acknowledgement. The authors would like to thank a Ph.D. student, Waixiang 
Cao for her assistance in the numerical experiments. 
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